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The topological quantum number Q of a superconducting or chiral insulating wire counts the 
number of stable bound states at the end points. We determine Q from the matrix r of reflection 
amplitudes from one of the ends, generalizing the known result in the absence of time-reversal and 
chiral symmetry to all five topologically nontrivial symmetry classes. The formula takes the form of 
the determinant, Pfaffian, or matrix signature of r, depending on whether r is a real matrix, a real 
antisymmetric matrix, or a Hermitian matrix. We apply this formula to calculate the topological 
quantum number of A'^ coupled dimerized polymer chains, including the effects of disorder in the 
hopping constants. The scattering theory relates a topological phase transition to a conductance 
peak, of quantized height and with a universal (symmetry class independent) line shape. Two peaks 
which merge are annihilated in the superconducting symmetry classes, while they reinforce each 
other in the chiral symmetry classes. 
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I. INTRODUCTION 

The bulk-boundary correspondence in the quantum 
Hall effect equates the number Q of occupied Landau lev- 
els in the two-dimensional bulk to the number of propa- 
gating states at the edge, which is the quantity measured 
in electrical conduction.— Thouless et al. identified Q as 
a topological quantum number;^ determined by an invari- 
ant integral of the Hamiltonian H{k) over the Brillouin 
zone. 

One-dimensional wire geometries can also be classified 
by a topological quantum number, which then counts 
the number of stable ("topologically protected") bound 
states at the end points. Examples exist in chiral insu- 
lators (such as a dimerized polyacetylene chain4) and in 
superconductors (such as a chiral p-wave wire^). In the 
former case the end states are half-integer charged soli- 
tons, in the latter case they are charge-neutral Majorana 
fermions. 

Following the line of thought from the quantum Hall 
effect, one might ask whether the number Q of these end 
states can be related to a transport property (electrical 
conduction for the insulators and thermal conduction for 
the superconductors). The basis for such a relationship 
would be an alternative formula for Q, not in terms of 
H{k)r-'— but in terms of the scattering matrix S of the 
wire, connected at the two ends to electron reservoirs. 

This analysis was recently carried out for the supercon- 
ducting p-wave wire,"i-^ which represents one of the five 
symmetry classes with a topologically nontrivial phase in 
a wire geometryj^^iH In this paper we extend the scat- 
tering theory of the topological quantum number to the 
other four symmetry classes, including the polyacetylene 
chain as an application. 

The outline is as follows. In the next section we 
show how to construct a topological invariant Q from 
the reflection matrix r (which is a subblock of S). De- 
pending on the presence or absence of particle-hole sym- 



metry, time-reversal symmetry, spin-rotation symmetry, 
and chiral (or sublattice) symmetry, this relation takes 
the form of a determinant, Pfaffian, or matrix signature 
(being the number of negative eigenvalues), see Table [B 
In Sec, mil we demonstrate that this Q indeed counts the 
number of topologically protected end states. The con- 
nection to electrical or thermal conduction is made in 
Sec. IIVI where we contrast the effect of disorder on the 
conductance in the superconducting and chiral insulat- 
ing symmetry classes. We conclude in Sec. |V]with the 
application to polyacetylene. 



II. TOPOLOGICAL QUANTUM NUMBER 
FROM REFLECTION MATRIX 

The classification of topological phases is commonly 
given in terms of the Hamiltonian of a closed systemi^ 
For the open systems considered here, the scattering ma- 
trix provides a more natural starting point. In an iV- 
mode wire the scattering matrix S* is a 2N x 2N unitary 
matrix, relating incoming to outgoing modes. The pres- 
ence or absence, at the Fermi energy Ep, of particle-hole 
symmetry, time-reversal symmetry, spin-rotation sym- 
metry, and chiral (or sublattice) symmetry restricts S 
to one of ten subspaces of the unitary group U{2N). In 
a one-dimensional wire geometry, five of these Altland- 
Zirnbauer symmetry classesAS. can be in a topological 
phase, distinguished by an integer-valued quantum num- 
ber Q. 

The symmetries of the scattering matrix in the five 
topological symmetry classes are summarized in Table [D 
For each class we have chosen a basis for the incoming and 
outgoing modes at the Fermi level in which the symmetry 
relations have a simple form. (In the next section we 
will be more specific about the choice of basis.) Notice 
that the chiral symmetry operation is the combination 
of particle-hole and time-reversal symmetry (if both are 
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TABLE I: Classification of the symmetries of the unitary scattering matrix S at the Fermi level in an A^'-mode wire geometry, 
and relation between the topological quantum number Q and the reflection submatrix r. For Z2 topological phases Q is given 
in terms of the sign of the determinant (Det) or PfafBan (Pf) of r. For Z topological phases the relation is in terms of the 
number 1/ of negative eigenvalues of r. 



present). 

Topological phases are characterized by a resonance 
at the Fermi level, signaling the presence of one or more 
quasi-bound states at the end-points of the wire with van- 
ishingly small excitation energy. (If the wire is supercon- 
ducting, these excitations are Majorana fermions.^) It is 
therefore natural to seek a relation between the topolog- 
ical quantum number Q and the reflection matrix, which 
is an iV X iV submatrix relating incoming and reflected 
modes from one end of the wire, 

-(:::)■ 

The wire has two ends, so there are two reflection ma- 



trices r and r' . Unitarity ensures that the Hermitian 
matrix products rr^ and r'r'^ have the same set of re- 
flection eigenvalues tanh^ A„ S (0,1), numbered by the 
mode index n — 1,2, .. .N. The real number A„ is the so- 
called Lyapunov exponent. The transmission eigenvalues 
T„ — 1 — tanh^ A„ = 1/ cosh^ A„ determine the conduc- 
tance G oc J2n wire. (Depending on the system, 
this can be a thermal or an electrical conductance.) The 
topological phases have an excitation gap, so the r„'s 
are exponentially small in general, except when the gap 
closes at a transition between two topological phases. A 
topological phase transition can therefore be identified 
by a sign change of a Lyapunov exponent.—"— 

The Lyapunov exponents are the radial variables of the 
polar decomposition of the scattering matrix, given by^S 
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in class D, 



in class Dili, 



in class BDI, 



in class AIII, 



Ql 
Ql 



in class CII, 



(2.2a) 
(2.2b) 
(2.2c) 
(2. 2d) 
(2.2e) 



in terms of a real diagonal matrix A — diag (Ai, A2, . . .) 
and complex unitary matrices Up (satisfying Up^ — U^), 
real orthogonal matrices Op (satisfying = 0| = 

Oj), and quaternion symplectic matrices Qp (satisfying 
Qp^ = Ql = ^yQp^y)- The matrices J^i crt ® ai ® 



• • - ©(Ti are block diagonal in terms of 2 x 2 Pauli matrices 
ai (with ao the 2x2 unit matrix). There are TV distinct 
A„'s in classes D, BDI, and AIII, but only N/2 in classes 
Dili and CII (because of a twofold Kramers degeneracy 
of the transmission eigenvalues). 

The transmission eigenvalues only determine the Lya- 
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punov exponents up to a sign. To fix the sign, we demand 
in class D and Dili that DetOp = 1, so Op £ SO(iV). 
Then the A„'s can be ordered uniquely as^^ |Ai| < A2 < 
A3 < • • • , so there can be at most a single negative Lya- 
punov exponent. In the other three classes there is no 
sign ambiguity since tanh A„ is an eigenvalue of the re- 
flection matrix r itself — which is a Hermitian matrix in 
classes BDI, AIII, and CII. There is then no constraint 
on the number of negative Lyapunov exponentSf^ 

If we start from an initial state with all A„'s positive, 
then the number Q of (distinct) negative Lyapunov ex- 
ponents in a final state counts the number of topological 
phase transitions that separate initial and final states. 
In class D this produces the relation Q = sign Det r from 
Refs. [TllfTsl . relating topological quantum number and 
determinant of reflection matrix. 

In class Dili the determinant of r is always positive, 
but we can use the Pfaffian of the antisymmetric reflec- 
tion matrix to count the number of negative A„'s, so we 
take Q = signPf r. [In view of the identity Pf XYX^ = 
(DetX)(Pfr), one has Pfr = {Bet Oi)Pf (A ^ icjy) = 
n„tanh A„.] 

In classes BDI and AIII the matrix signature Q = i'{r) 
of the Hermitian matrix r gives the number of negative 
eigenvalues, equal to the number of negative A„'s. In 
class CII we take Q — ^^{r) to obtain the number of 
distinct negative A„'s, because each eigenvalue is twofold 
degenerate. 

These topological quantum numbers are defined rel- 
ative to a particular reference state, chosen to have all 
positive Lyapunov exponents. We would like to relate 
Q to the number of end states at zero excitation energy, 
and then chose a reference state such that this relation- 
ship takes a simple form. This is worked out in the next 
section, with the resulting expressions for Q given in Ta- 
ble H 



III. NUMBER OF END STATES FROM 
TOPOLOGICAL QUANTUM NUMBER 

We consider first the superconducting symmetry 
classes D and Dili and then the chiral symmetry classes 
BDI, AIII, and CII. The symmetry class D was treated 
in detail in Ref. [ll| and is included here for completeness 
and for comparison with class DHL 



A. Superconducting symmetry classes 

Electron-hole symmetry in a superconductor relates 
the energy- dependent creation and annihilation opera- 
tors by ^\E) = ^{—E). Since therefore 7^ = 7 at i? = 0, 
an excitation at zero energy is a Majorana fermion, equal 
to its own antiparticle. The end states in symmetry 
classes D and Dili are socalled Majorana bound statesi^ 
In the open systems considered here, where the supercon- 
ducting wire is connected to semi-infinite normal-metal 



FIG. 1; Superconducting wire (S) connected to a normal- 
metal lead (N) which is closed at one end. A bound state at 
the Fermi level can form at the NS interface, characterized 
by a unit eigenvalue of the product tn^ns of two matrices of 
reflection amplitudes (indicated schematically by arrows). 



leads, the end states are actually only quasi-bound states, 
but they still manifest themselves as a resonance in a con- 
duction experiment.— 

The topological quantum number in class D should 
give the parity of the number Af of Majorana bound 
states at one end of the wire: Af is even {Q = 1) in the 
topologically trivial phase, while JV is odd (Q — —1) in 
the topologically nontrivial phase. In class Dili all states 
are twofold Kramers degenerate so Af is to be replaced 
by Af/2. 

Let us now verify that the determinant and Pfaffian 
expressions for the topological charge in Table H] indeed 
give this bound state parity. We transform the quasi- 
bound states into true bound states by terminating the 
normal-metal lead at some distance far from the normal- 
superconductor (NS) interface (see Fig.[T|). For the same 
purpose we assume that the superconducting wire is suf- 
ficiently long that transmission of quasiparticles from one 
end to the other can be neglected. The reflection matrix 
tns from the NS interface is then an TV x iV unitary ma- 
trix. The number of modes N = 2M is even, because 
there is an equal number of electron and hole modes. 

The condition for a bound state at the Fermi level is 



Det (1 - tnt-ns) = 0, 



(3.1) 



where tn is the reflection matrix from the terminated 
normal-metal lead. In the electron-hole basis the matrix 
tn has the block-diagonal form 



C/n 

, m 



(3.2) 



The matrix J7n is an M x M unitary matrix of electron 
reflection amplitudes. The corresponding matrix for hole 
reflections is because of particle-hole symmetry at the 
Fermi level. 

The reflection matrix from the NS interface has also 
off-diagonal blocks, 




(3.3) 



Particle-hole symmetry relates the complex reflection 
matrices The ~ r*f^ (from electron to hole and from hole 



to electron) and r-g, 
from hole to hole). 



' hh 



(from electron to electron and 
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Class D 



A unitary transformation, 



(3.4) 



produces a real reflection matrix r — r*. This is the 
socalled Majorana basis used for class D in Table HI 
The determinant is unchanged by the change of basis, 
Det tns = Det r. 

The condition (13.11) for a bound state reads, in terms 
of r, 



Det (1 + O^r) = 0, 



(3.5) 



with On = — r^^NS^^ an orthogonal matrix. The number 
J\f of bound states is the number of eigenvalues —1 of 
the 2M X 2M orthogonal matrix Oj^r, while the other 
2M — TV eigenvalues are either equal to +1 or come in 
conjugate pairs e^*'^. Hence Det On?" — (—l)"^ and since 
Det On = 1 we conclude that Detr — (—1)"^, so indeed 
the determinant of the reflection matrix gives the bound 
state parity in class D. 



2. Class Din 

Time-reversal symmetry in class Dili requires 

Ans = i^yrm = -^nS' (3-6) 



with E 
now define 



y 



ay. Instead of Eq. p.4|) we 



(The matrix T,y acts on the spin degree of freedom, hence 
it commutes with J7, which acts on the electron- hole de- 
gree of freedom.) In this basis r = r* is still real, as re- 
quired by particle-hole symmetry, while the time-reversal 
symmetry requirement reads r = —r"^. This is the basis 
used for class Dili in Table HI 

The Pfaffians in the two bases are related by Pf r = 
(Detrj)(Pf Ans) = (-l)^/*Pf Ans- (Each electron and 
each hole mode has a twofold Kramers degeneracy, so the 
total number of modes N is an integer multiple of four.) 
The relation can be written equivalently as 



Piir = Pf A 



NS- 



(3.8) 



This identity is at the origin of the factor i appearing 
in the class Dili expression for the topological quantum 
number in Table HI 

The condition (13.11) for a bound state can be rewritten 



(3.9) 



Det (An -r) = [Pf (An - r)]^ = 0, 



where An = ft{i'Syr]^)H'^ , as well as r, are antisymmet- 
ric orthogonal matrices. In App. |X] we show that the 



multiplicity TV of the number of solutions to Eq. p.9p 
satisfies 

(_l)A^/2 ^ (Pf AN)(Pfr). (3.10) 

Since, in view of Eq. p.2p . 

Pf Aat = (Detf7)|Pf (iI]j,C/t)p = (-1)^/^, (3.11) 

we conclude that (— l)^/'*Pfr = Pfir gives the parity 
of the number TV/2 of Kramers degenerate bound states. 
This is the topological quantum number for class Dili 
listed in Table HI 

B. Chiral symmetry classes 

In the chiral symmetry classes BDI, AIII, and CII we 
wish to relate the number iy{r) of negative eigenvalues 
of the refiection matrix r to the number of quasi-bound 
states at the end of the wire. As before, we transform 
these end states into true bound states by terminating 
the wire and assume that the transmission probability 
through the wire is negligibly small (so r is unitary). 
While in the superconducting symmetry classes we could 
choose a normal metal lead as a unique termination, in 
the chiral classes there is more arbitrariness in the choice 
of the unitary refiection matrix ro of the termination. 

Since refiection matrices in the chiral classes are Her- 
mitian (see Table H| , we can decompose 



UoSuoUq, 



1 



JV-rio 







-1„ 



(3.12) 



(3.7) where C/q is an TV x TV unitary matrix, no — i^iro), and 



1„Q is an no x uq unit matrix. (Unitarity restricts the 
eigenvalues to ±1.) Similarly, 



(3.13) 



with iy{r) = ni. 

Time-reversal symmetry with (without) spin-rotation 
symmetry restricts the unitary matrices Uq and Ui to the 
orthogonal (symplectic) subgroup, but to determine the 
number of bound states we only need the unitarity. 

The condition Det (1— ror) — for a zero-energy bound 
state takes the form 



Det(l-5„„C/5„,C/^) =0, 



(3.14) 



with U = U^Ui. We seek the minimal multiplicity TV 
of the solutions of this equation, for arbitrary U. (There 
may be more solutions for a special choice of U, but these 
do not play a role in the characterization of the topolog- 
ical phase.) 

We divide U into four rectangular subblocks. 



U 



N—na,N—ni 



M 



N—riQ ,ni 



M, 



no.N—ni 



(3.15) 
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where M„ m is a matrix of dimensions n x m. Since 



in view of unitarity of U, the bound state equation p.l4p 
simphfies to 



Det 







Mn- 




(3.17) 



For any matrix Mn^m with n < m there exist at 
least m — n independent vectors v of rank m such that 
Mn^mV = 0. Therefore Eq. p.l7p has at least \no+ni—N\ 
independent solutions, hence 



JV^ \i^{r) + i^{ra) ~ N\. 



(3.18) 



This is the required relation between the topological 
quantum number Q — i^(r) (in class BDI, AIII) or 
Q = ^i^{r) (in class CII) and the minimal number of 
bound states Af at one end of the wire, for arbitrary ter- 
mination of the wire. In the special case of termination 
by a reflection matrix tq = — Itv =^ ^^(''o) — N, the re- 
lation takes the simple form A/" = Q (in class BDI, AIII) 
and A/'^ 2Q (in class CII). 

So far we considered one of the two ends of the wire, 
with reflection matrix r. The other end has reflection 
matrix r' = —r [see Eq. (|2.2p ]. so J^(r') = N — i/(r). Ter- 
mination of that end by a reflection matrix Tq produces 
a minimal number A/"' of bound states given by 



M' = Hr) - Hr'o)\- 



(3.19) 



For Tq = Ijv ^ >^{r'o) — we have the simple rela- 
tion A/"' = Q (in class BDI, AIII) and M' = 2Q (in class 
CII). The (minimal) number of bound states at the two 
ends is then the same, but in general it may be differ- 
ent, depending on how the wire is terminatedj^i^ This 
arbitrariness in the chiral symmetry classes is again in 
contrast to the superconducting classes, where Majorana 
bound states come in pairs at opposite ends of the wire. 



IV. SUPERCONDUCTING VERSUS CHIRAL 
SYMMETRY CLASSES 

As a first application of our general considerations, we 
contrast the effect of disorder and intermode scattering 
on topological phase transitions in the superconducting 
and chiral symmetry classes. We focus on the symme- 
try classes D and BDI, which in the single-mode case 
are identical, so that the effect of intermode scattering is 
most apparent. 

In both these classes there is particle-hole symmetry, 
which implies that we can choose a basis such that the 
Hermitian Hamiltonian H satisfies 



We assume for simplicity that the N right-moving and 
left-moving modes all have the same Fermi velocity vp- 
To first order in momentum p = —ihd/dx the Hamilto- 
nian then takes the form 



+ Uo[iAix) (g) + iB{x) (g) 



^C{x)®a,X (4.2) 



with Iat the N x N unit matrix. The N x N matrices A 
and B are real antisymmetric, while C is real symmetric. 
(For A^ = 1 this model Hamiltonian was used in Ref. ll.) 

The Hamiltonian (|4.2[) respects all the symmetries 
present in class D, but in class BDI the additional chiral 
symmetry requires 



UxHux = —H. 



(4.3) 



This implies that the matrix B = in class BDI. 

The transfer matrix M. relates the wave function "^{x) 
at the two ends of the disordered wire (of length L): 
vI'(L) = A^\1>(0). At the Fermi level (zero energy) M 
follows upon integration of the wave equation H"^ — 
from X = Q to X ^ L, 



A4 = Texp|^^ j da; ^-AqIat (8) CTj; 
+ Uo[A{x) (To + iB{x) ®(jy- C{x) ® Ox 



(4.4) 



The symbol 7" indicates the ordering of the noncommut- 
ing matrices in order of decreasing x. 

The Pauli matrices in Eq. (|4.4p define a 2 x 2 block 
structure for the 2A^ x 2A^ transfer matrix. The N x N 
reflection matrix r and transmission matrix t follow from 
this block structure by solving 



(4.5) 



The reflection matrix gives the topological quantum num- 
ber, Q — sign Det r in class D and Q = i'{r) in class BDI. 
The transmission matrix gives the conductance 




G = GoTrtt^ 



(4.6) 



H* 



H. 



(4.1) 



In class D this is a thermal conductance (with Go = 
TT^k'^To/Qh, at temperature tq), while in class BDI this 
is an electrical conductance (with Go = 2e^//i). 

We model a disordered wire in class D by taking a 
Gaussian distribution (zero average, unit variance) of the 
independent matrix elements oi A(x) , B (x) , C (x) , piece- 
wise constant over a series of segments of length SL <^ L. 
In class BDI we use the same model with B = 0. 

In Fig.[2]we plot the conductance and topological quan- 
tum number as a function of Ao for different values of C/q, 
calculated in class D and BDI for a single realization of 
the disorder. A change in Q is accompanied by a peak 
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FIG. 2: Conductance G (top panels) and topological quantum number Q (bottom panels) in the superconducting class D (left 
panels) and the chiral class BDI (right panels). The black and blue curves are calculated from the Hamiltonian (|4.2|l . for a 
single disorder realization in a wire with N = 5 modes. The red dotted curve show the universal line shape (|4.7p of an isolated 
conductance peak. Energies Aq and Uo are measured in units of hvp/SL for SL = L/10. 



in G, quantized at Gq if the topological phase transi- 
tions are well separatedJ^ The difference between the Z2 
superconducting topological phase and the Z chiral topo- 
logical phase becomes evident when conductance peaks 
merge: In the superconducting class D the conductance 
peaks annihilate, while in the chiral class BDI a maxi- 
mum of N conductance peaks can reinforce each other. 

Also shown in Fig. [2] is that a single isolated conduc- 
tance peak at Aq = has the same line shape as a 
function of (5 = (Aq - Ac)/r, 

Go 



Gpcak(^) 



cosh^ S ' 



(4.7) 



in both the superconducting and chiral symmetry classes. 
(The width T of the peak is not universal.) We have 
checked that the Une shape in the other three symmetry 
classes also has the same form (|4.7p . so this is a general 
statement. One cannot, therefore, distinguish the Z2 and 
Z topological phases by studying a single phase transi- 
tion. This is a manifestation of the super-universality of 
Gruzberg, Read, and VishveshwaraJ^ 



V. APPLICATION TO DIMERIZED POLYMER 
CHAINS 

We conclude with an application in a physical sys- 
tem. Such an application was given for the supercon- 
ducting symmetry class D in Ref. fill , so here we con- 
centrate on the chiral classes. We consider a dimerized 
polymer chain such as polyacetylene, with alternating 
long and short bonds, described by the Su-SchriefFer- 
Heeger Hamiltonian.^^ This is a tight-binding Hamilto- 
nian, which in the continuum limit takes the form of the 



class BDI Hamiltonian (|4.2[) .^^ Our goal is to obtain the 
Z topological quantum number of N coupled polymer 
chains from the reflection matrix. 

The single-chain electronic Hamiltonian is2^"— 



Nl 



n=l 



tn+l,n — to — Ct{Un+l — Un) , 



(5.1a) 

(5.1b) 



with to and a nearest-neighbor (real) hopping constants 
and c„ the electron annihilation operator at site n. (The 
spin degree of freedom plays no role and is omitted.) Chi- 
ral (or sublattice) symmetry means that H 1-^ —H ii 
Cn I— > — c„ on all even-numbered or on all odd-numbered 
sites. We take Nl even, so that the chain contains an 
equal number of sites on each sublattice. 

Following Jackiw and SemenofF'^® we ignore the atomic 
dynamics, assuming that the electrons hop in a pre- 
scribed atomic displacement field of the dimerized form 
Un — (— l)"uo+<5u„. Disorder is accounted for by random 
displacements Sun, chosen independently on N parallel 
chains. Nearest neighbors on adjacent chains are coupled 
by an interchain hopping constant iintcr, which we take 
non-fluctuating for simplicity. 

The reflection and transmission matrices r and t were 
computed from the Hamiltonian (15. ip via the transfer 
matrix, as outlined in App.|B] In Fig.[3]wc show the topo- 
logical quantum number Q (equal to the number j/(r) of 
negative eigenvalues of the Hermitian reflection matrix 
r), as well as the electrical conductance G — GoTrtt^ 
(with Go = 2e^/h). These two quantities are plotted as 
a function of the dimerization parameter mq, to illustrate 
the topological phase transition, but unlike the excitation 
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FIG. 3: Conductance (black dotted line, left axis) and topo- 
logical quantum number (blue solid line, right axis) of = 3 
coupled polymer chains (each containing A''^ = 300 sites). 
These curves are calculated from the reflection and transmis- 
sion matrices, obtained from the Ifamiltonian (|5.ip for to = 1, 
a = \, and tintcr = 0.1, for a single realization of the random 
5it„'s (having a Gaussian distribution with zero average and 
standard deviation 0.2). The red dotted curve shows the uni- 
versal line shape (|4.7|) of an isolated conductance peak. 



gap Ao in a superconducting wire this is not an externally 
controllable parameter. 

The case iV = 3 plotted in Fig. [3] is a Z topological 
phase, and each change in the topological quantum num- 
ber is accompanied by a peak of quantized conductance. 
The lineshape again has the universal form (|4.7p . 
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Appendix A: Calculation of the number of end 
states in class Dili 



We wish to prove that the multiplicity M of the number 
of solutions of the bound state equation (|3.9p satisfies Eq. 
p.lOp . for arbitrary antisymmetric orthogonal matrices 
and r of dimension N x N, with N = 2M and M an 
even integer. 

We use that any antisymmetric orthogonal matrix can 



be factorized as = On«S„Om 



terms of orthogonal matrices On and Ons- These fac- 
torizations relate a Pfaffian to a determinant, Pf An ~ 
DetON, Pir = DetONS- 

We seek the multiplicity JV of the number of solutions 

of 

[Pf (An - r)f = ^ [Pf (iEy - OiY^yO^)] ^ = 0, (Al) 

with O = On^ns an orthogonal matrix. 

We consider the secular equation for the twofold de- 



generate eigenvalues z„ of the matrix iT,yOiT,yO'^ , 
= Det (z - iEj^OiSyO'^) = Det (ziE^^ + OiT^yO'^) 



[Pf {ziJ:y + OiT.yO^)Y 



M 



= Pf (ziJ^y + OiT.yO'^) = cY[iz- 2„) = 0. (A2) 



The value c = 1 of the prefactor follows by sending z to 
infinity. By filling in z = we find that 



M 



Pf (OiSyO'^) = DetO = Jl 



(A3) 



n=l 



The 7V/2 bound state solutions have z„ — —1, the 
remaining M — J\f /2 solutions have either z„ = 1 or con- 
jugate pairs z„ = e'^^'^. Hence 



M 



n ^« - (-1)-^/' = Det O = (Pf A^.)(Pf r), (A4) 
as we set out to prove. 



Appendix B: Calculation of the topological quantum 
number of a dimerized polymer chain 

To simplify the notation we outline the calculation of 
the topological quantum number for the case A'^ = 1 of 
a single polymer chain, when the transmission matrix r 
is a scalar and we may take Q = ^(1 — Q') with Q' = 
signr G {—1, !}• (The multi-chain case, with Q = z/(r) G 
{0, 1,2,... A^}, is analogous.) 

From the tight-binding Hamiltonian (|5.1|) we directly 
read off the zero-energy transfer matrix A4 in the site 
basis, 



Mn 




(Bl) 
(B2) 



The normalization factors in Eq. (|B1I) have been inserted 
so that the current operator has the site-independent 
form X = Uy. 

To obtain the scattering matrix we need to transform 
from the site basis to a basis of left-movers and right- 
movers, in which the current operator equals rather 
than ay. This change of basis is realized by the matrix 
f2 from Eq. 



(B3) 
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Multiplying the transfer matrices we find for the entire 
chain (containing an even number of sites Nj^): 



M^Mn,Mn,-i---M2Mi= 1^ , (B4) 



2X 1x2-1 + 1 



with the definition 



Nl/2 



Nl/2 ] \ ^2n+l,2n 
^2n,2n-l 



(B5) 



(B6) 



We obtain the reflection and transmission amplitudes 
from M with the help of Eq. (|4.5p . The result is 



(B7) 



I- 2X 
r=- t 



so the topological quantum number is given by 
Q' = sign(l-X2) 

fNL/2 Nl/2 \ 

= Sign n ^2n,2n-l " H ^2n+l,2n • (B8) 



ri=l 



If all hopping constants are close to to > we may 
simplify this expression to 



'Nl/2 



Q' = sign ^ [t2«,2n-l - t2«+l,2n] 



(B9) 



In the absence of disorder, when t2n,2n-i = ~ 2auo, 
i2ri+i,2n = ^0 + 2auoj this reduces further to Q' — 
— sign auo , so we recover the original criterion that the 
dimerized polymer chain has bound states at the ends if 
the weaker bond is at the endi^ 
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